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Auto CPAP 

CROSS REFERENCE TO RELATED APPLICATIONS 

The present application claims the benefit of Swedish patent 
5 application No. SE 0003531-1 which was filed on October 2, 2000 and 
is incorporated herein by reference. 

FIELD OF THE INVENTION 

10 The present invention relates to method and an apparatus for the 
detection and treatment of disordered breathing during sleep, in 
particular to a method and apparatus employing an artificial neural 
network. 

15 BACKGROUND OF THE INVENTION 



|; = U.S. Patent No. 5,953,713 (Behbehani et al.), incorporated herein by 

rU reference, discloses a method for treating sleep disordered breathing 

comprising measuring a respiration-related variable at an interface 
2 0 placed over a patient's airway coupled to a pressurized gas, feeding 
cepstrum data obtained from the respiration related variable(s) into 
an artificial neural network trained to recognize patterns 
characterizing sleep disordered breathing; supplying pressurized gas 
to the patients airway in response to recognition of the artificial 
25 neural network of sleep disordered breathing. The sampling frequency 
of the pressure transducer's output disclosed in the preferred 
embodiment is 512 Hz. A Fourier transform is calculated every 1/16 
second using a 32 sample values window. 
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Another aspect of frequency analysis is that, on the one hand, the 
precision is proportional to the number of input data but that, on the 



other hand, the response time is correspondingly increased. While 
high precision is welcome since rather small changes in breath 
pattern can be detected, a slower response increases the risk of 
progressive deterioration of the airway aperture, and thereby more 
severe respiratory disturbance before the patient is aroused. Other 
methods of detecting sleep disorder are based on breath-by-breath 
analysis 

Alternatively, if adequate treatment is not installed, the patient will 
be aroused in a more extended time perspective. 

OBJECTS OF THE INVENTION 

It is an object of the present invention to provide an improved 
method for automatically supplying continuous positive airways 
pressure to a patient 

It is another object of the present invention to provide an automatic 
continuous positive airways pressure apparatus (ACPAP) which lacks 
at least some of the drawbacks of state-of-the-art apparatus. 

Additional objects of the invention are evident from the following 
short description of the invention, the attached drawings illustrating a 
preferred embodiment, the detailed description thereof, and the 
appended claims. 

SUMMARY OF THE INVENTION 

The present invention is based on the insight that a direct analysis of 
the flow signal is more specific than an analysis of disordered breath, 
in particular flow limited breath, based on frequency analysis. 




According to the present invention is provided an automatic 
continuous positive airways pressure apparatus (ACPAP) in which the 
air flow from a CPAP or other system providing positive air pressure 
to a patient is measured for calculation of a number of parameters 
specific to the signal. The set of parameters comprises cepstrum 
coefficients and energy content, and is selected to indicate an apneic 
event of breathing during sleep, such as apnea, hypoapnea, and flow 
limitation. Data for these parameters collected from a large number 
of patients were used to train an artificial neural network to teach the 
system the variation ranges of the parameters for subsets of patients 
under a number of circumstances. The result from the artificial neural 
network is obtained as a low-dimensional grid of nodes in which each 
respiration type is represented by trajectory or a subsets of nodes. 
A trajectory for a normal breath looks very different from that of a 
breath during disturbed sleep. 

If breathes symptomatic of a condition of disturbed sleep are 
detected the CPAP pressure is increased. In contrast, CPAP pressure 
is reduced in a normal condition. 

Thus, according to the present invention is disclosed a method for the 
detection and treatment of disordered breathing during sleep 
employing an artificial neural network in which data related to 
25 breathing gas flow are analyzed in an artificial neural network. 

Specifically, a method according to an embodiment of the present 
invention comprises the following steps: 

placing a mask with a tube over a patient's airway, the mask being 
30 in communication with a source of a pressurized breathing gas 
controlled by a CPAP, thereby establishing a respiratory circuit; 
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periodically sampling the gas flow in the circuit; 
periodically calculating values for one or several parameters 
distinctive of a breathing pattern; 

periodically feeding the parameter values to an artificial neural 
5 network trained to recognize breathing patterns characteristic of 
sleep disordered breathing; 

analyzing the parameter values in the neural network; 
controlling pressurized breathing gas supply in response to the 
output from the neural network. 

10 

It is preferred to feed said parameter values to the network at a 

a 

i.n frequency of from 2 Hz to 30 Hz, preferably of about 20 Hz. It is 

preferred for said parameters to comprise cepstrum coefficients and 
5 energy slope. 

CO 15 

According to a first preferred aspect of the invention the artificial 
neural network is trained with data collected from a large number of 
W patients. The data will have been collected from patients differing in 

□ many aspects: sex, age, body weight, breath pattern, etc. In 

20 addition, variants of sleep disordered breathing such as those 

occurring preferentially in the back position, those occurring during 
particular stages of sleep, and those occurring under the influence of 
drugs or alcohol need to be addressed. 



25 Such data are advantageously collected in sleep laboratories in which 
the state of sleep is followed as well as the type and severity of the 
breathing disturbance is monitored by use of a polysomnography 
system. The collected data form a primary database. During the 
training of the artificial neural network the data is quantified under 

30 formation of a small secondary dedicated database which can be 
stored in a ACPAP. Thus, according to the present invention, a 



dedicated secondary database obtained from a primary database 
comprising data collected from a large number of persons is stored in 
the ACPAP. 

5 According to a second aspect of the invention it is preferred to 
periodically sample the gas flow during breathing. 

The ANN comprises a number of nodes representing sets of training 
data. Each note reflects a state or an incident (feature). Neighboring 
nodes represent incidents of small geometric distance. In the same 
way as in training an incident vector is extracted for each flow data 
sample. The Euclidean distance from the incident vector to each node 
is calculated. The node in closest proximity to the vector is associated 
with it. Sequences of incident vectors are followed as sequences of 
nodes in the artificial neural network. It can be said that a sequence 
of nodes is the response of the network. Thus a trajectory in the 
geometric structure of the network (response) is followed rather than 
in the parameter space. The fact that the dimension of the network 
most often is smaller than the parameter space is of advantage since 
calculation thereby is simplified. The response from the network 
forms the basis for distinguishing between apnea, hypoapnea and a 
normal breathing state and thus, for CPAP pressure control. 

The invention thus is based on the use of an artificial neural network 
25 (ANN) of Kohonen-map type (associative memory; T. Kohonen, Self- 
Organization and Associative Memory, 2 nd Ed., Springer Verl., Berlin 
1987) for detecting apnea or apnea-like episodes. The ANN is trained 
with data obtained from a number of patients in a sleep laboratory. 
The readily trained ANN forms a global (universal) structure of data 
30 stored in a non-volatile memory in an ACPAP. In use the breathing 
pattern of a patient forms trajectories (traces) in the ANN. A normal 
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breathing cycle forms a closed trajectory. A trajectory deviating 
from normal is indicative of a breath disturbance. The ANN is 
structured in way so as to make certain areas represent initial stages 
of apnea. The passage of a trajectory through such an area or 
5 touching its border indicates that the amount of air provided to the 
patient should be increased so as to re-establish normal breathing. 
Once breathing has been normalized the adduced amount of air is 
reduced to normal, i.e., to the pre-established base value. 

The artificial neural network is trained in two phases described in P. 
Brauer and P. Knagenhjelm, Infrastructure in Kohonen Maps, Proc. 
IEEE ICASSP, Glasgow 1989. 

The purpose of the analysis is to extract, from the series of air flow 
rate measurements, values of the parameters chosen to classify and 
detect apneic and hypoapneic states. In each single analysis the 
parameters are made to form an incident or feature vector on which 
all training and decision-making is based. All sample values are 
individually analyzed in preparation for a quick response to changes 
in flow which are typical forewarnings of an apneic or hypoapneic 
state. 

According to a third preferred aspect of the invention linear predictive 
coding is used to analyze the parameter values fed to the neural 
25 network. A linear predictive coding analysis comprising four 

parameters is carried out for all samples. In particular, the so-called 
A-parameters from the analysis are converted to cepstrum 
parameters for optimal correlation between parameter distance and 
conceptual distance, that is, so-called associativity. 
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According to a fourth preferred aspect of the invention the 
prediction error in calculating linear predictive coding is used as a 
basis for the parameter next in line. The error is filtered to counteract 
short-term variations and normalized with the total energy of the 
5 analytical window. 

For calculations of energy a larger window than for the linear 
predictive coding analysis is used. The energy of the latest windows 
can be used to calculate a line the inclination which describes a trend. 
10 The difference in trend is used as a further parameter. Thus, 

according to a fifth preferred aspect of the invention, the inclination 
of a trend line calculated from measurements and is used as a 
parameter. 

15 According to the present invention is also disclosed an apparatus for 
the detection and treatment of disordered breathing during sleep for 
use with a CPAP, the apparatus including a probe for sampling 
breathing air flow data, in particular on inhalation, and an artificial 
neural network for analyzing, directly or indirectly, said data to 

20 control breathing air pressure. 

According to the present invention is also disclosed a CAPAP 
comprising a probe for sampling breathing air flow data, in particular 
on inhalation, and an artificial neural network for analyzing, directly 
25 or indirectly, said data to control breathing air pressure. 

Further variations of the present invention are disclosed in the 
following detailed description of a preferred embodiment thereof 
illustrated in a drawing. 
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DESCRIPTION OF THE DRAWING 

The invention is illustrated by a drawing comprising several figures, 
showing: 

5 

Fig. 1 depicts a block diagram regarding parameter extraction; 
Fig. 2 depicts a feature map response in regard of example 3; 
Fig. 3 depicts a feature vector; 

Fig. 4 depicts a feature map response in regard of example 4. 

10 

DESCRIPTION OF A PREFERRED EMBODIMENT 

|ij EXAMPLE 1. General 

M 15 PRIMARY DATA ANALYSIS. The purpose of the analysis is to extract, 
is from the series of air flow rate measurements, values of the 

3 parameters chosen to classify and detect apneic and hypoapneic 

\*\ states. In each single analysis the parameters are made to form an 

Q incident or feature vector on which all training and decision-making is 

20 based. All sample values are individually analyzed in preparation for a 

quick response to changes in flow which are typical forewarnings of 

an apneic or hypoapneic state. 

Incident vector parameters . LPC-Cepstrum. For each sample a Linear 
2 5 Predictive Coding (LPC) four- para meter analysis is carried out. The 
so-called A-parameters from the analysis are converted to cepstrum 
parameters for optimal correlation between parameter distance and 
conceptual distance, that is, associativity. The term cepstrum 
introduced by Bogert et al. in connection with echo time series 
30 analysis designates the inverse Fourier transform of the logarithm of 
the power spectrum of a signal. The transformation of a signal into its 




cepstrum is a homo-morphic transformation, see A.V. Oppenheim 
and R.W. Schafer, Discrete-Time Signal Processing, Prentice Hall, 
Englewood Cliffs, NJ, 1989. Residual. The error of prediction in 
calculating LPC is used as a basis for the following parameter. The 
5 error is filtered to oppose short-term variations, and is normalized 
with the total energy for the analytical window. Energy slope. For 
calculations of energy larger windows are used than for LPC analysis. 
The energy at the most recent windows is used to calculate a line the 
slope of which describes a trend. Difference in trend. The difference in 
10 trend is used as a further parameter. 

PARAMETERS. To detect an apneic event (i.e. central/obstructive 
apnea, hypoapnea, and flow limitations) a model is used to 
5 characterize typical qualities and features of the flow-signal during 

:ij the event. The parameters of the model is chosen with the aim to be 

!°° 15 as distinct, unambiguous, and informative as possible. The set of 

parameters shall respond to typical apneic events that are readily 

2 detected by physicians. In addition to be sensitive to apneic events, 

[ y , it is important that the parameters shall be insensitive to features 

p irrelevant to the task of detecting apneic events. 

20 THE FEATURE VECTOR. The values of the parameters are compiled to 
form a vector, below named the Feature Vector. Each time the flow- 
signal is measured (sampled), the values of the Feature Vector are 
extracted. This means that if the flow signal is measured f s times per 

second, and N parameter values are needed in the model, the data 
25 rate is increased from /, to N- f s samples per second. 

Prior to the extraction of parameter values, the flow signal is 
differentiated (high-pass filtered) to avoid the influence of the mean 
signal value. The mean will vary with patients and/or hardware and 
do not contribute in the classification of apneic events, and is 



• 10 • 



therefore removed. Each N -dimensional Feature Vector can be 
regarded as one point in a -dimensional signal space. 

TRAINING THE NETWORK. An Artificial Neural Network (ANN) is 
iteratively trained to organize groups or clusters of Feature Vectors 
5 with similar properties. The self organizing process known as 
Kohonen's Self-Organizing Feature Map [1-2] has shown great 
capability of performing this task. 

The number of clusters is defined prior to the training and is 
10 determined by the required resolution of the ANN. The training is 
initiated by a set of M clusters, randomly positioned in the N - 
dimensional signal space. 

The database used for training is formed by compiling the Feature 
15 Vectors from a large number of patients with various sleep disorders 
and at all stages of sleep. During the training, each input Feature 
Vector is compared to each cluster to find the one with best 
resemblance to the input vector. This cluster is voted winner, and is 
adjusted towards the input vector. In addition, all other clusters 
20 within a neighborhood to the winner in another domain, the so-called 
map-space are adjusted towards the input vector. The map-space is 
usually of low dimension containing one node for each cluster in the 
signal-space. The nodes are arranged in hexagonal or a square 
lattice, and the Euclidian distance between them defines their internal 
25 relation. A node's neighborhood is usually defined by a neighborhood 
function and contains the set all nodes in the beginning of the training 
whereas only a few (or none) are considered neighbors at the end. 
The further away a node is to the winner in the map-space, the less 
the corresponding cluster in the signal-space is adjusted towards the 
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input vector. Thus all adjustments are done in the signal space, 
while the rules of adjustments are defined in the map-space. 

The training time is predetermined, and an annealing function is used 
5 to "freeze" down the system causing only small adjustments at the 
end of the training. The neighborhood function creates correlation 
between the signal-space distance and the map-space distance 
allowing classification to be performed in the (low dimensional) map- 
space, rather than in the more complicated signal-space. The method 
10 described above is known as "unsupervised learning", i.e. there is no 
need to use classified data in the training procedure described above. 
Classification of data into various apneic events is a tedious task. 

When the ANN is readily trained, the clusters will represent 
15 ^features of the input flow signal including normal breathing, 

hypoapnea, flow-limitations, and apnea (provided these features are 
represented in the database used for training). The response of the 
ANN is proportional to the signal distance between the input signal 
and all the clusters. See figure 2. Often this output is of less interest 
20 in the case of classification. The output is instead used to find the 
node with best resemblance to a classified input, such as normal 
breathing and apneic events. This is known as the labeling phase in 
the design of the ANN. Classified Feature Vectors are presented for 
the ANN, the output is observed, and the node giving the highest 
25 output is labeled with the presented feature. The actual output 
thereafter is the label rather than the response value. 

The set of clusters are now stored in the memory of the APAP to be 
used in runtime mode. Patient flow-data is analyzed exactly the same 
30 way as done in the training phase to extract the values of the 

parameters used in the model i.e. the Feature Vector. The vector is 



then presented to the network that will produce the output label 
(classification) which is used by the flow-control logic. 

EXAMPLE 2. Algorithm 

DATA ACQUISITION. 

5 Let the flow signal be a digitized version of the analog flow-signal 
sampled at f s samples/second, giving a sequence of samples 

x„ * = 0,1,... , 

where is short for x(i-T) i.e. the sample at time instant i-T and 

g T=\if s . 

CH io PREPROCESSING. To reduce the influence of individual patient 
5 variations and to facilitate classification stability, the signal should 

12 pass a device to remove the signal mean. Any kind of steep edge 

high pass filter can be employed, thus the ideal differentiator is used 

l3 

t g for simplicity. The output from the differentiator, d, (and the input to 

\j 15 the parameter extractor) will then be 

'ssa? 

M- rf, = x, - x M where / = o,l,... , and *_,=0 

PARAMETERS. The cepstrum coefficients have shown to well model 
the frequency content of the signal using only a few parameters (low 
order model). In addition, the dynamics of the cepstrum coefficients 
20 facilitate quantization of the parameters. Often the parameters are 
weighted to produce parameters with similar variances. The cepstrum 
coefficients are derivatives of the so called A-polynomial calculated by 
standard Linear Predictive Coding (LPC). 



25 



As the cepstrum coefficients used do not hold information about the 
signal energy, the cepstrum will be augmented with a parameter to 
account for the long term (say 10 seconds) energy variations. This 



parameter shall be insensitive to the absolute level of the flow signal 
and only reflect the relative fluctuations. 

Cepstrum. T o calculate p model parameters at time k , the last Wz 
5 input samples are used (windowing). The sample, d k , is predicted to 

be 
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d k =~X a / ' d k-i ■ 



/=i 



Thus the prediction error signal is 

e k =d k -d k =-^a l -d k _ { 



10 The task is to find the set, a, which minimizes the energy of the 

prediction error signal (i.e. finding the values of a = [a,, to make 

d k as similar to d k as possible) over all samples within the window. 

The optimal solution is found solving the matrix equation 

R a = -r , 

15 known as the Yule-Walker or normal equation, using for instance the 
Cholesky algorithm. 

The vector a constitutes the so called A-polynomiai which is 
transformed into a set of cepstrum coefficients c using the following 
algorithm: 



c 0 =0 



The cepstrum coefficients c X9 ...,c P are used as the p first coefficients 
in the Feature Vector. 



Parameter of Energy Trend (PET) . The energy from a number of 
windows are used to calculate a trend for the energy values. If the 
trend indicates increasing energy levels, PET is set to zero. If the 
trend is decreasing, the point where the trend is crossing the time- 
axis is calculated. A non-linear transform of this value form the PET 
coefficient. HBK figure 

The PET is added to the Feature Vector. The complete Feature Vector 
is thus composed according to Figure 3. 

FEATURE MAP GEOMETRY AND DEFINITIONS. Let the Mmap nodes 
be denoted 

m /5 / = 0 ? ...,M -1 . 

Most often the nodes are arranged in a square (2-dimensional) grid. 
The distance between two map nodes / and j , is denoted d u and 

defined as the squared Euclidian distance (Z 2 norm) between them 

=£ 2 (m,,m y ). 

Let the input Feature Vector representing sample x k be denoted y k . 

The map response in node / for feature k , s ik , is defined as: 
■i « 

S lk =exp£ d*/(P + \)} f 
where the signal space distance df k is defined as 

l=\ 

and w i is some suitable weight function. 

ANNEALING FUNCTION. The task of the annealing function is to 
obtain an equilibrium at the end of the training. The principle is that 
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large adjustments are allowed in the beginning of the training 
whereas only small (or zero) adjustments are allowed at the end. 
How the decrease incorporated is not critical. Linear, exponential, and 
even pulsating [4] decay schedules are proposed in the literature. 

5 

INITIALIZATION. Traditionally, all data driven clustering schemes, 
including ANNs, employ random positioning of the clusters in the 
signal space, by assigning (small) random numbers to the 
parameters. The actual values are not important as long as they are 
10 not identical. The ordering of the clusters is also at random. 

TRAINING. The iterative algorithm adjust all clusters after each input 
Feature Vector, y k / presented. The direction of the adjustment is 

towards y kr and how much is determined partly by the annealing 

15 function, partly by the neighborhood function. The adjustment 
formulae for cluster y k at time instant r + l is: 

y * (' + 1) =y * (0 + Y k - (0(z - y k (0) , 

where 

Y k '(t)=f(0'g k (t) 

20 and f(t) is the annealing function and g(t) is the neighborhood 
function. Various suitable functions are discussed in [3]. 



EXAMPLE 3. Hard decision A 

25 Let 64 map nodes be arranged in an 8x8 square grid and numbered 0 
to 63 from the top left to the low right corner. Thus for example the 
map distance £> 01 = 1, D 02 =4 and d 09 = 2 ■ 
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A large database is recorded containing flow- measures from several 
patients during all phases of sleep. The recordings are performed at 
20 Hz and stored on a memory disk. The database will contain normal 
sleep breathing, flow limitations, snoring, yawning, coughing, various 
apneic events, but also mask leakage and other artifacts. 

The database is analyzed sample for sample. The 20Hz flow-signal is 
first passed through an ideal differentiator. A rectangular window of 
180 samples is used to form basis for extracting 4 cepstrum 
coefficients (c 1? ...,c 4 ) and the PET parameter. Thus the Feature Vector 

is a 5-dimensional vector with values extracted every 50ms. 

Samples are collected from the database in a random manner as long 
as the training proceed. The number of iterations, T, is determined by 
the size of the database, but as a rule of thumb, 10-30 iterations per 
sample may be an adequate number. 

The Euclidian distance to all clusters y / / = 0,...,M-1 are calculated 

and the cluster closest to the Feature Vector is voted winner and is 
denoted y, . 

In calculating the distance, the following weight function is employed: 



The neighborhood function will allow all clusters to be adjusted at all 
times (i.e. the size of neighborhood is not decreased in time), but will 
penalize clusters far away from the winner y } . 

g(t) = e- 2D » ,W 

The annealing function,/^), follows a linear decay schedule 

/(/)=0.2. a -//T) * = 0,...,T 



w 



= [1 2 3 4 l] 
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After the training phase, the map is presented with known features 
such as normal breathing, flow limitation signals etc. For each event 
the response, S, in each node of the map is calculated (see, Fig. 2). 



Regions with high response for normal breathing are labeled as 
normal regions; regions reacting for flow-limitations are labeled as 
flow-limitation area and so forth. 

It is then decided which nodes that shall represent events needed to 
respond to; i.e. an alarm signal is passed to the pressure control 
system. In Fig. 2, the map response to an input corresponding to flow 
limitations is depicted. In this case, nodes X,Y,Z will probably be 
labeled as a flow limitation region. 



This example is similar to the one described in example 1, but here a 
Hamming window of 180 samples is used instead of the rectangular. 
Furthermore the cepstrum coefficients are weighted to have 
approximately the same variances, whereas the PET parameter is 
given twice the variance of the cepstrum. This will give the PET 
parameter a little more importance than the rest of the parameters. 
The following weight function is employed: 



The training is carried out exactly ass before, and the map response 
for a typical flow-limitation signal is shown in Fig. 4. In this case, 
nodes X,Y,Z will probably be labeled as a flow limitation region. 



S ik =exp(-^/i>), 



EXAMPLE 4. Hard decision B 



= [l 1.41 1.5 4 2] 
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EXAMPLE 5. Soft decision 

In this example, the map nodes are given a number 0, 2, 5, or 10 as 
labels to indicate the seriousness of the classification result. Thus 
instead of presenting an alarm signal on/off to the pressure control 
5 system, the number is passed on. If for instance 0 reflects normal 
regions and 10 reflects apnea, the numbers can be integrated to form 
an overall breathing status classification. If the level is very high, 
rapid increases in pressure is allowed, low levels allow for a pressure 
decrease, and intermediate levels result in a slow increase of the 
io pressure. 

EXAMPLE 6. Down sampling 

If the capacity of the processor do not allow for all calculations 
described above, the LPC calculation can be decimated by a factor 
two or four. The number of samples within the analysis window must 
15 then be reduced so that the time span of the window is still about two 
breathing cycles. The resolution of the map response will not suffer 
from this. 

EXAMPLE 7. Pressure regulation and runtime mode 

The pressure control system will increase the pressure one step of 
20 0.125 mbar if the ANN response is positive for 50 samples in one 
sequence. The pressure will decrease with one step of 0.125 mbar if 
the ANN response is negative for 300 samples in one sequence. The 
pressure will not be changed if the ANN response is changed during 
those sequences. 

25 In a runtime mode, a sample of the flow signal is analyzed as 

described above (i.e. extracting the Feature Vector), and presented to 
the map, now stored in a memory bank in the APAP-unit. There is no 
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need to calculate the exponent in the expression for map response, 
as the function is monotonic. 



EXAMPLE 8. Program Code 

5 // function ann detect 

// COPYRIGHT 7C) 2000 PePe Research 
// GOTHEBORGH, SWEDEN. 

#include <stdio.h> 
10 #include "ann_detect.h" 
#include "apne_func.h" 

int ann_detect(const float new_sample) { 

15 float corr[ 10]; 

int ana_cfim = 5; 

r3 int win = 180; 

5 float apol[5]; 

:= float alpha; 

~ 2 0 float dvec[5]; 

'■11 int respx, respy; 

'jjj int resp_en, resp_ann; 

i!y calc_corr(new_sample, win, ana_dim+l, corr); 

M 2 5 

/* Check levels, and slope */ 
resp_en = reg_r0(corr[0], &dvec[4]); 



dvec 
dvec 
dvec 
35 dvec 



lpc(corr, ana_dim, apol, &alpha); 
a2cep(apol, ana_dim, dvec); 



*= 1.41; 
*= 1.5; 
*= 4" 



map_resp(dvec, ana_dim, &respx, &respy); 
resp_ann = apne_dec(respx, respy); 

40 if ( resp_ann == 1 || resp_en = = 1 ) 

return 1; 

return 0; 



45 // functions for apne detection 

f // COPYRIGHT (C) 2000 PePe Research 
// GOTHEBORGH, SWEDEN. 

50 #define MAP_DIM_X 8 

#define MAP_DIM_Y 8 

#define MAP_SIZE 64 

#define DATA DIM 5 

#define NHEAD 256 
55 #define NMB_SAVE 100 

#define MaxRO 180 



40 
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#define MaxRO_ 2 90 

#define SLOPE 0.006 
#include <stdio.h> 
#include <math.h> 
#include "map.h" 
#include "annpar.h" 



float distance ( const float *x, const float *y, int d ) 

10 int i; 

float dist, t; 



^t x i o tt£ [0]; 

15 for ( i = l- i<d; i++ ) { 
t = x[i]-y[i]; 
dist += t*t; 



20 } 



return dist; 



i= ; void calc_corr( float const rec_sample, int const win, int const 

ana_dim, float* corr) 

! S 25 < 

w static float samp_win[180]; 

static int first = 0; 

m static float prev_sample=0; 

lO float diff; 

u 30 int i,j; 

~ /* calc diff of new and prev sample */ 

'■■i diff = rec_sample - prev_sample; 

'»H prev_sample = rec_sample; 

ru 3 5 

hJ /* Insert new sample first, move all the other */ 

r=t for ( i=win-l; i > 0; i— ) { 

Cl samp_win[i] = samp_win[i-l]; 

samp_win[0] = diff; 

for ( j = 0; j < ana_dim; j++){ 
corrlj] = 0; 

for (i = 0; i < win-j; i++K 

corr[j] += samp_win[i]*samp_win[i+j]; 

> } 

50 /* Add an offset to avoid near /0 operations */ 
corr[0] +=1; 



55 int reg_r0(float cur_r0) 

^ static float rf[MaxR0]; 
static int first = 1; 
static int index = 0; 
60 int i; 

static float prev_r0 = 0; 



float k,m, rfO, dec; 

rfO = 0.9*prev_r0 + 0.1*cur_r0; 
prev_rO = rfO; 

k = 0.5*(rf0 - rf[MaxRO-l]): 

m = (rfO + rf[MaxR0/2] + rf[MaxR0-l])/3.; 

/*update the saved old values */ 
10 for ( i=MaxR0-l; i > 0; i-- ) { 
rf[i] = rf[i-l]; 
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rf[0] = rfO; 



/* second rule check the regression line */ 
if (k > -1) 

dec = 1000; 
else 

20 dec = -MaxR0_2*((m/k)+l); 

Q *out = exp(-SLOPE*dec); 

! /* first rule, if energy less then 15 make alarm */ 
K 2 5 if (cur_r0 < 15 ){ 

return 1; 

S > 

i;fl /* second rule, if regression line crosses zero line too soon make 

u 30 alarm */ 
:. if ( dec < 50 ){ 

return 1; 

} 

:H return 0; 

W 35 } 

\\ 

a void lpc(float* const acf, int const ana_dim, float* aCoeffs, float* 

i e ± alpha) 
{ 

4 0 float err; 
int i, j; 

float refls[5]; 

float sum, refl, tmp; 

4 5 err = acf[0J; 

for (i = 0; i < ana_dim; i++) { 
if (err <= 0) { 
refls[il = 0; 
aCoeffs[i] = 0; 
50 } else { 

sum = 0; 

for (j = 0; j < i; j++) 

sum += aCoeffs[j] * acf[i - j]; 
refl = (acf[i + 1] - sum) / err; 
55 refls[i] = refl; 

forfi = " ■ - 



(J = 0; i < i/ 2; j++) { 
tmp = aCoeffs[jl; 
aCoeffsM] -= refl * aCoe 



iCoeffs[i - j - 1]; 
aCoeffs[i - j - l]~-= refl * tmp; 

60 } 

if (i & 1) 



en 
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aCoeffs[i / 2] -= refl * aCoeffs[i / 2]; 
aCoeffsR] = refl; 
err *= 1 - refl * refl; 

for (i = 0; i < ana_dim- i++) { 
aCoeffs[i] = -aCoeffs[i]; 

10 ^alpha = sqrt(err); 
> 

void a2cep(float* const a, int const ana_dim, float* cep) 

15 { 

int I, j, ind; 
float sum; 

cep[0] = -l*a[0]; 
20 for (I = 2; I < ana_dim+l; I++) { 
sum = l*a[l-l]; 
for (j=2; i < = I; j++) { 
ind = l-j+1; 

sum = sum + a[j-2]*cep[ind-l]*ind; 

25 } 

cep[l-l] = -sum/I; 



m 
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void map_resp(float* const dvec, int const ana_dim, int* respx, int* 
respy) 
{ 

int xm, ym, ind, i; 
35 float dist; 

float min_dist = 1000000; 
for ( xm = 0; xm < MAP_DIM_X: xm ++ ) { 
for ( ym = 0; ym < MAP_DIM_Y; ym ++ ) { 
ind = (xm + (ym*MAP_DIM_X))*DATA DIM; 
4 0 if ( (dist = distance(dvec, &map[ind], DATA_DIM)) < 

min_dist ){ 

*respx = xm; 
* respy = ym; 
min_dist = dist; 



45 



50 



> } 



int apne_dec( int const respx, int const respy, float* d_out) 

* static float mean_x = MAP_DIM_X/2; 
static float mean_y = MAP_DIM_Y/2; 
55 static int first = 0; 

static float resp_hist[NMB_SAVE]; 
static int n_hist_x = 0; 
static int n hist_y = 0; 
float rad_fiTt; 
60 float rad inst; 

double dx, dy, dfx, dfy; 



10 
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int current_in_area = 0; 
int new_in = 0, i; 

dx - (respx-area_x_mid)*(respx-area_x_mid); 
dy = (respy-area_y_mid)*(respy-area_y_mid); 
rad_Jnst = sqrt(dx + dy); 

mean_x = lp_a*mean_x + ( l-lp_a)*respx; 
mean_y = lp_a*mean_y + (l-lp_a)*respy; 

dfx = fmean_x-x_cen)*(mean_x-x_cen); 
dfy = (mean_y-y_cen)*(mean_y-y_cen); 
rad_filt = sqrt(dfx +dry); 

15 /*check if the current response is in the area */ 
if ( area_rad > radjnst ) new_in = 1; 

for ( i = NMB SAVE-1; i > 0; i~ ) { 
resp_hist[T] = resp_hist[i-l]; 

20 } 

resp_hist[0] = new_in; 

for M = 0; i < nb resp hist; i++ ) { 

if ( 1 == resp3ist[iT) current_in_area++; 

25 > 

if ( mean_x > det_x_min && 
mean_y > det_y_min && 
mean_x < det_x_max && 
mean„y < det_y_max) 
30 /* standard case, when detecting a simple square area of the 

filtred response */ 
return 0; 
else if (rad > rad_filt ) 

/* within circle calc, otherwise same as above */ 
35 return 0; 

else if (resp_in_area < current_in_area ) 

/*instant resp totaly over the given number */ 
return 0; 
else 

4 0 return 1; 

}; 



4 5 EXAMPLE 9. Patient treatment 

Equipment . Breas CPAP pvlO prototype (Breas AB, Molndal, Sweden) 
with internal flow measurement; PSG system, EMBLA™ 
Polysomnography (Embla hf, Reykjavik, Iceland); PC with artificial 
50 neural network software for sleep disturbance detection. 

Patients . 6 males aged from 31 to 60 years, suffering from sleep 
disorders of various kind. 




Measurement. The patient was set up with all sensors from the PSG 
system: EEG, EOG, EMG, oxygen saturation, pulse, nasal flow 
thermistor measurement, body position. With this setup it is possible 
to determine sleep stages, arousals, sleep apnea, sleep hypoapnea, 
5 and other sleep related events. The CPAP provides air through a tube 
via a nasal mask to the patient. Via the CPAP communication 
interface the PC was connected to the CPAP. The PC software can 
read air flow values from the CPAP and set new pressure set points 
on the CPAP by the communication interface. The information was 

10 exchanged at a rate of about 20 Hz. The PC program feeds the flow 
values into the artificial neural network (ANN). The output from the 
ANN is entered into a pressure regulation algorithm (PRA). The 
pressure regulation algorithm calculates a new pressure set point and 
activates the new value in the CPAP. The output from the artificial 

15 neural network and the pressure setpoint is read by the PSG system. 
Evaluation . In a normal CPAP titration (sleep disorder analysis) a 
CPAP is connected to a patient during sleep. The CPAP pressure is 
adjusted during the night so as to put the patient in a state with no 
indications of sleep disorders. This pressure is the one used in the 

2 0 CPAP treatment. The patient's need for various CPAP pressures can 
be seen with the PSG system breath-by-breath. The required CPAP 
pressure varies depending on sleep stage, body position, etc.. The 
data from the auto CPAP test was evaluated in the PSG system in the 
same manner as for CPAP titration by a physician used to evaluate 

25 patients receiving CPAP treatment. Thereby the correlation between 
the detection of sleep disorder by the ANN and the analysis in the 
PSG system could be determined. The good correlation obtained 
indicated that the AFN reacted correctly. 
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